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ABSTRACT 

We report the discovery of a six- month-long mid-infrared transient, SDWFS- 
MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO 
Deep Wide-Field Survey Bootes field. The transient, located in a 2; = 0.19 low 
luminosity (M[4.5] ~ —18.6 mag, L/L^ ~ 0.01) metal-poor (12 + log(0/H) 2^ 
7.8) irregular galaxy, peaked at a mid- infrared absolute magnitude of M[4 5] ~ 
—24.2 in the 4.5yum Spitzer/IRAC band and emitted a total energy of at least 
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]^q5i gj^gg x]2g optical emission was likely fainter than the mid-infrared, although 
our constraints on the optical emission are poor because the transient peaked 
when the source was "behind" the Sun. The Spitzer data are consistent with 
emission by a modified black body with a temperature of ~ 1350 K. We rule 
out a number of scenarios for the origin of the transient such as a Galactic star, 
AGN activity, GRB, tidal disruption of a star by a black hole and gravitational 
lensing. The most plausible scenario is a supernova exploding inside a massive, 
optically thick circumstellar medium, composed of multiple shells of previously 
ejected material. If the proposed scenario is correct, then a significant fraction 
(~ 10%) of the most luminous supernova may be self-enshrouded by dust not only 
before but also after the supernova occurs. The spectral energy distribution of 
the progenitor of such a supernova would be a slightly cooler version of rj Carina, 
peaking at 20-30/im. 

Subject headings: galaxies: irregular — infrared: galaxies — supernovae: general 
— supernovae: individual (SDWFS-MT-1, SN 2007va) 



Introduction 



The first truly large variability survey of extragalactic sources at mi d-infrared (naid-IR , 
3.6-8.0/xm) wavelengths is the Spitzer peey, Wide-Field Surv ey (SDWFS: lAshby et al.ll2009[ ) 
of the NOAO Deep Wide-Field Survey (jjannuzi fc Dey 1119991 ) Bootes field. The survey spans 
the years 2004-2008 with fou r epochs, covers 8 deg^, and contains variability statistics for 
nearly half a million sources ( iKozlowski et al.ll2010al ). While the majority (~ 76%) of the 
extragalactic variable objects in the mid-IR are active galactic nuclei (AGNs), there is room 
for serendipitous discoveries such as supernovae (SNe). 

In general, SDWFS was not expected to be interesting for SNe searches. First, the 
contrast between SNe and their hosts is relatively poor in the mid-IR, essentially because of 
differences in the effective temperatures. If we match the bolometric luminosities of an SN 
and a typical host galaxy spectral energy distribution (SED), the SN contributes ~ 70% of 
the V-band flux but only ~ 8% of the 3.6/im flux. Second, with only four epochs of data, we 
expected SDWFS would add little to our understanding SNs even if one were detected. Third, 
the expected rates are very low. Typica l Type la (Hp) SN p eak at M[36] ~ ~ —18.4 
(—18.1) mag (e.g.. iKrisciunas et al.ll2004j : iMaguire et al.ll2010l ). so SDWFS can detect them 
in a single epoch only for redshifts 2: ^ .07 (0.06). Given the SN e la rates from the S DSS- 
II Supernova Survey (IDilday et al.ll2010[ ) and SNe Type II rate (IHoriuchi et al.ll2009[ ). we 
would then expect to detect only of order one Type la SN and of order seven Type II SNe 
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in the SDWFS survey even if we could safely search for them at the 5a detection level 
for a single epoch ([3.6] = 19.1 mag). In practice, controlling the false positive rates for 
identifying variable sources means that the variability selectio n criteria must be signi ficantly 
more conservative than a 5a peak at the detection limit (see iKozlowski et al.ll2010al ). so we 
did not expect to detect any SNe in SDWFS. 

There is a subclass of Type I In SN e that are far more luminous, reaching M[3 g] ~ ~ 
—20 mag (e.g., iMannucci et al.ll2003l ). which we could detect up to 2; ~ 0.15. The leading 
theory for these SNe uses collisions between two n iassive shells of materi a l to efficiently 
radiate the kinetic energy of the SN in the optical (jSmith fc McCrayl 120071 . IWoosley et al. 
20071 ). The fast moving SN ejecta collides with an outer, slower moving shell on scales of order 
10^ AU where the optical depth is high enough to produce a well-defined photosphere but 
low enough for the thermal energy from the shock to be radiated before adiabatic expansion 
converts it back into kinetic e nergy. But these SNe are rare, representing only < 1/160 
of the local Type II SNe rate (IMiller et al.l 120091 ). which more than balances the increased 
detection volume. 



Here we investigate the nature the brightest mid-IR transient in SDWFS. Over a 6 
month period, this z ~ 0.2 source radiated ~ 10^^ ergs at a black body temperature of 
~ 1350 K. In Section |2] we present the available UV, optical and mid-IR data along with 
a Keck spectrum of the host galaxy. Then, in Section |3] we consider a range of possible 
scenarios for producing this source and conclude that a simple variant of the models for 
the hyper-luminous Type Iln SNe is the best explanation for this source. The paper is 
summarized in Section HI Throughout this paper we use a standard ACDM model with 
{QA,QM,^k) = (0.73,0.27,0.0) and h = i/o/100 = 0.710. All magnitudes are in the Vega 
system. 



2. Data 



2.1. Spitzer/IRAC Data 



SDWFS-MT-1 (aka SN 2007va; IKozlowski et al.ll2010bl ). where MT means mid-IR tran- 
sient as a parallel notation to using OT for optical transient, is the most significantly 
variable source in the entire SDWFS field. It corresponds to the SDWFS source SDWFS 
J142623. 24+353529.1. The [3.6] and [4.5] light curves were strongly correlated (r = 1) 



^To derive basic parameters of the transient and its host galaxy, we used the Cosmology Calculator 
(|Wrightll2006[ ): http : //www . astro . ucla . ed.u/~wr ight/ CosmoCalc . html 
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and the variability amplitude was 76 standard deviations from the mea n for its average 
magnit ude. A detailed definition of these variability criteria is given in iKozlowski et al. 
( l2010a( ). The SDWFS survey consists of four epochs, taken 3.5 years, 6 months, and 1 



month apart. The object is not detected i n any of the four IRAC bands in the first epoch 
(also known as the IRAC Shallow Survey, lEisenhardt et al.l 120041 ) on 2004 January 10-14 
(JD' = JD-2450000 ^ 3016.5). The SDWFS single epoch 3cr detection limits of [3.6] = 19.67 
mag, [4.5] = 18.73 mag, [5.8] = 16.33 mag, and [8.0] = 15.67 mag (4 arcsec apertures cor- 
rected to 24 arcsec; Vega mag), set our upper limits on the object brightness in the first 
epoch. The observed part of the transient peaked in the second epoch taken 2007 August 
8-13 (JD' ^ 4322.5) at [3.6] = 15.93 mag and [4.5] = 15.61 mag (Figures [U and [2]). It 
was slightly fainter 6 months later in epoch 3 (2008 February 2-6, JD' ^ 4500.5), with 
[3.6] = 16.09 mag and [4.5] = 15.76 mag. Then, in a matter of a month, the object faded by 
almost 2 magnitudes to [3.6] = 17.99 mag and [4.5] = 17.72 mag in epoch 4 on 2008 March 
6-10 (JD' ^ 4533.5). We present all the available Spitzer data, in Tabled! 



IRAC flFazio et al.l 120041 ) simultaneously observes the [3.6]/[5.8] or [4.5]/[8.0] bands, 
and in SDWFS the [4.5]/[8.0] observations were taken 40 seconds after the [3.6]/[5.8] ones. 
Each epoch consists of three 30 second exposures in each band, taken no more than 2 days 
apart, and the images for the whole mosaics were taken within 5 days. We confirmed that 
the transient was present in the individual artifact-corrected basic calibrated data (CBCD) 
frames corresponding to the position of the transient. The mosaic image for each epoch is 
then a combination of three to four of these fram es, and there is no indication of problems 
in the SDWFS reductions 



see 



Ashbv et al.ll2009h 



2.2. Keck Spectrum 



We obtained a Keck/LRIS ( lOke et al.lll995[ ) spectrum of the transient's host galaxy on 
2010 March 12, long after the transient was gone. The spectrum (shown in Figure |3]) shows 
clear emission lines of Ha, H/3, H7, H5, He, and HC, [OH] at 3727A, [O HI] at 4363A, 4959A 
and 5007A, [Nelll] at 3869A and 3968A. The redshift of the host galaxy is z = 0.1907, so 
the luminosity distance is 920 Mpc, the distance modulus is 39.82 mag, and the angular 
scale is 3.15 kpc/arcsec. Before measuring the line fiuxes, we corrected the host g alaxy 
spectrum for Galactic extinction (E(B— V) = 0.014 mag) using ISchlegel et al.l (119981 ). We 
do not detect the [Nil] 6584A line, which implies log([NII]/Ha) < —1.5. Combining it with 
log([QIII] 5007A/H/3) = 0.78, we conclud e that the host is an irregular HII galaxy (see, e.g., 
Terlevich et"aDll99ll : iKniazev et aDl2004h . The lack of [Nil] and [NeV] 3426A fines and the 
clear detection of the 4363A line rule out the presence of an AGN. 
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We followed the prescription of iPilyugin fc ThuanI ( 120071 ) to derive the gas-phase oxygen 
abundance of the host. We measured the fluxes of the [Oil] 3727A, [OIII] 4363A, [OIII] 
4959A, [OIII] 5007A and H/3 486lA emission lines to find a metaUicity of 12 + log(0/H) = 
7.85±0.13. We assumed an electron density of rie = 100 cm~^, but the value does not change 
significantly with changes in the electron density (e.g., 12 + log(0/H) = 7.87 if we use = 
10 00 cm~^ inst e ad). N ext, we verified this metaUicity measurement using the prescriptions 
of llzotov et al.l (l2006l ) to get 12 + log(0/H) = 7.79. As a cross check, we also used the 
oxygen-to- hydrogen flux rati o R23 = ([Oil] + [0III])/H/3 = (F^.^y^ + FA4959 + i^A5007)/i^A486i 
method of iPagel et al.l ( 119791 ) and the lPettini fc Pagell ( 120041 ) 03N2 method to derive oxygen 
abundances. The clear detection of the 4363A line and log([NII]/Ha) < —1.5 strongly 
suggests that we should use th e "l ower" metaUicity branch in the R23 method. We use the 
relations from ISkillmanI (119891 ) and lYin et al.l ( l2007l ) with our measured log(-R23) = 0.94, to 
get oxygen abundances of 12 -|- log(0/H) = 7.78 and 7.80, respectively. The 03N2 method 
gives 12 + log(0/H ) < 8.0. All these es timates are approximately 10% of the solar value 
(assuming 8.8 from iDelahaye et al.ll2010[ ). 



2.3. UV and Optical Data 



In order to understand the origin of this mid-IR transient, we searched for any available 
UV and optical data for this area of the sky in the relevant time span. 

The host galaxy w as detected in the original NOAO Deep Wide-Field Survey (NDWFS, 
Jannuzi fc Il999[ ) with = 23.98 ± 0.05 mag {XFx = 6.64 x 10^^ ergs/s), R = 
22.78 ± 0.05 mag (1.10 x 10^^ ergs/s), and I = 22.58 ± 0.09 mag (8.39 x 10^^ ergs/s). 
It was unresolved in 1.2 arcsec (FWHM) seein g, which puts an upp er limit on the size 
of the galaxy of 3.8 kpc. In archival GALEX (iMorrissey et al.l l2007l ) data, the host has 
NUV= 24.30 ±0.17 mag and no detection in FUV. The GALEX source is formally 0'.'87 from 
the SDWFS/NDWFS source, but there are no other blue sources within 20 arcsec. Both 
the colors of the host and its size are consistent with a blue, low lumi nosity, low rnetalli city 
galaxy. We used the template model of an irregular galaxy from lAssef et al.l (I2OIOI ) to 
estimate the mid-IR magnitudes of the host. Fitting the GALEX NUV, and NDWFS Bw, 
R, and /-band magnitudes of the host galaxy, and assuming z = 0.1907, the expected 
apparent (absolute) IRAC magnitudes are [3.6] = 20.87 (-18.56) mag, [4.5] = 20.69 (-18.57) 
mag, [5.8] = 20.52 (-20.02) mag, and [8.0] = 18.37 (-21.45) mag. This is consistent with 
the failure to detect the host in the first SDWFS epoch. We estimate rest frame host 
luminos ities of Mr = —16 .11 mag and Mr = —16.75 mag, making this an L/L^ ~ 0.01 
galaxy (IKarachentsevI l2005l ) . 
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Although the field was monitored by at least 5 optical variability surveys, the peak of 
the transient occurred when the field was "behind" the Sun, so almost all the optical data 
corresponds to the periods before epoch 2 and after epoch 4. The Bootes field was observed 
seven times between JD' = 3800 and 4620 with the 8kx8k Mosaic CCD camera on the 2.4 
meter MDM Hiltner Telescope. The two most interesting observations took place 90 days 
before the second and 84 days after the fourth SDWFS epoch, bracketing the period of the 
transient. In these 180 second exposures, we detected no fiux from the transient (or the 
host) at a 3cr level of R = 21.32 mag (4.2 x 10^^ ergs/s) and 21.15 mag (4.9 x 10^^ ergs/s), 
respectively. 



The Cat alina Sky Survey (CSS jLarson et al.ll2003l )/Catalina Real-Time Transient Sur- 
vey (CRTS, iDrake et al.l l2009l ) provides light curves for the period prior to epoch 2 and 
overlapping epochs 3 and 4. Each epoch of four 30 second exposures on the 0.7 meter 
Catalina Schmidt Telescope reaches an r-band magnitude limit of ~ 20 mag. These data 
provide the strongest limit on the optical-to-mid-IR fiux ratios, where a stack of 12 images (3 
epochs), gives a 3a upper limit on the transient of r > 20.93 mag, corresponding to 7.6 x 10^^ 
ergs/s (large dark green symbol to the right in Figure [2]), in the period between epochs 3 
and 4. 



The QUEST (IDjorgovski et al.ll2008l ) survey further limits the optical emission before 
epoch 2 and after epoch 4. There are 19 epochs each consisting of 1 to 4, 60-100 sec exposures 
with the 1.2m Samuel Oschin telescope in a wide red/IR filter that we calibrated to i-band. 
The most interesting data were taken ~ 40 days prior epoch 2, where we stacked six adjacent 
QUEST epochs (JD' = 4276.71-4286.76) to obtain 3a detection limit of i > 20.70 mag 
(7.7 X 10^^ ergs/s, large orange symbol in Figure [2]). 

Sir nilarly, by coadding t he high cadence RAPTOR (RAPid Telescopes for Optical Re- 



sponse, IWozniak et al.l 120091 ) observations of the field, each of which is a 30 sec exposure 
using the RAPTOR-P array of four 200-mm Canon telephoto lenses, we obtained 15 upper 
limits at various epochs between 2004 and 2008, before and after the transient (Figure [2]). 
Calibrating the unfiltered observations to an i?-band equivalent using Tycho 2 stars, we find 
no evidence for the transient, with Cr > 18 mag (3o"). 



Only the northern station of the All Sky Automated Survey (ASAS) data (jPojmahski 



19971 ) obtained data during the peak of the transient. Combining 7 high quality ASAS 
images during the transient (JD' = 4322.75-4525.11), we obtain a weak upper limit on the 
transient magnitude of V > 16.3 mag {3a), corresponding to 6.2 x 10'^'^ ergs/s. 
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Discussion of the Transient 



The steeply faUing, roughly power law SED of the transient (Figure H]) is suggestive of 
thermal emission by a source cooler than normal stars. If we fit the SED as dust emission 
using an emissivity of Q oc A~^, we find good fits with /3 ^ 1 a s is typical of the mid-IR 



emission by warm dust (e.g. iSeki fc YamamotdllQSOl : lDraindll98ll : iGraham fc Meiklelll986l ). 
The temperature of T 1350 K varies little between epochs (Tab l e [J) and is below typical 



dust destruction temperatures of T > 1600 K (e.g., iDwekl 1 19831 : iGuhathakurta fc Draine 
19891 ). Table [1] summarizes these estimates of the temperatures and luminosities. A zeroth 
order estimate (ignoring the dust emissivity) of the radius R of the dust photosphere is 

R = 8000 ( J ^V'Vlgijiy AU. ,1) 

V1043 ergs/sy V ^ / 

We are assuming that the shell is optically thick and R corresponds to the radius where the 
(optical) optical depth is of order unity. If the dust is optically thin, as the case for typical 
dust echos, there is no characteristic radius or temperature since photons can be absorbed 
and emitted at any distance from the SN. We discuss these issues in more detail below. 
For the second and third epochs of SDWFS, with total luminosity L ^ Q ^ 10^^ ergs/s and 
temperature T 1350 K, we find the radius to be i? ~ 11000 AU (0.05 pc). The total 



energy radiated in the IR is then 



rad 



10 ergs. The weak optical limits mean that a 



comparable amount of energy might have been radiated in the optical. 



We used the more sophisticated radiation transport model DUSTY (jivezic fc Elitzur 



19971 ) to check this estimate. For incident black body spectra with temperatures 10,000- 
20,000 K, shell geometries with constant density, and dust temperatures at the inner pho- 
tosphere edge of 1300-2000 K, we could match the epoch 2 and 3 SEDs, including a weak 
optical limit, using shells dominated by graphite, a visual optical depth of r ~ 1-3, and a 
radius R ~ 10, 000-40, 000 AU. Pure silicate dusts were unable to fit the data well, produc- 
ing the wrong spectral shape for reasonable r and inner dust temperatures. These results 
are consistent with our simple estimates and strongly suggest that there was an optical 
counterpart, albeit with significant extinction [R > 21.7 mag). 

Before outlining our preferred explanation of the transient as a self-obscured supernova, 
we briefiy outline the hypotheses we rejected. 



3.1. Rejected Hypotheses 



First, the transient cannot be explained by a Galactic star. The arguments against are 
as follows: (1) based on the Keck spectrum, we know that the source is spatially coincident 
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with an irregular galaxy; (2) few stars can brighten by at least 5 magnitudes (by more than 
a hundred times) in the mid-IR and then stay at this level for 6 months before fading; and 
(3) most flaring stars do not flare in the mid-IR but in the UV. There are exceptions to (2) 
and (3), and these will generally be accreting systems where the disk can alter the spectral 
properties. However, the NDWFS field is a high latitude {b ~ 67°) field with negligible 
Galactic extinction and unlikely to contain young (protostars) or massive (Be stars) with 
these characteristics, particularly if they must also have no optical counterpart in quiescence 
despite needing to be nearby sources because of the short path length out of the disk. 

The second possibility is that we are really observing an AGN. First, we note that 
there are no signs of AGN activity in the host spectrum. The narrow lines in a galaxy 
are on relatively large scales, and should show line ratios averaged over decades (or more) 
of AGN activity, so the lack of a spectrum during the transient does not invalidate using 
them to argue against an AGN. Moreover, the mid-IR colors at the peak of the transient, 
[3.6] — [4.5] = 0.3 mag and [5.8] — [8.0] = 0.0 mag (epoch 2), are incon sistent with that of an 



AGN fe.g.. lStern et al.ll2005l : lAssef et al.ll2010HKozlowski et al.ll2010af ). While the composite 
SED of an AGN and its host can have these colors, we know from the host flux limits in the 
first epoch that the SED at the peak should be dominated by the AGN, and a pure AGN at 
2: ^ 0.2 should ha ve mid-IR colors closer to [3.6] — [4.5] ~ 0.8 mag and [5.8] — [8.0] ^1.1 mag 
( Assef et al. We also considered the p ossibility of a tidally d isrupted star in the context 

of the super-Eddington accretion models of iLoeb &: Ulmerl (119971 ). While these models have 
an inflated, relatively cool envelope around the black hole, the expected temperatures of 
T ^ 7000(Mbh/10^ Mq)1/^ K for a black hole of mass Mbh are still too high. A black hole 
with low enough mass (Mbh ~ 150 Mq) to match the temperature (T ^ 1350 K) would be 
unable to produce the observed luminosity. The Eddington luminosity for such a low mass 
is only LEdd ~ 2 x 10"^° ergs/s. 

The third possibility is that it is a 7-ray burst (GRB). We can rule it out as direct 
emission from a burst because th e time scale (6 m onths) is too long and the SED is wrong 



for direct beamed emission (e.g., ISari et al.lll999[ ). It cannot be indirect emission, where a 



dusty region absorbs and reradiates emission from the jet, because a relativistic jet expands 
past the scale of the photosphere (~ 10^ AU) too quickly (~ 60 rest-frame days) to produce 
the observed event duration. 



Finally, this cannot be gravitational lensing of a background object by the host galaxy 



see 



Wambsgansd l2006l for a review). The time scales are such that it would have to be 



microlensing by the stars in the lens galaxy, but there are no known mid-IR sources that 
are sufficiently compact to show such a large, and apparently achromatic, degree of magni- 
fication. Moreover, almost any potential source would have to show a still stronger degree 
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of optical magnific ation that would disagree even with our weak optical limits (see e.g., 
Morgan et aPboioi ). 



3.2. An Obscured Supernova Explosion 

Our working hypothesis is that S DWFS-MT-l is a supernova analogous to SN 2006gy, 
an extremely luminous Type I l n SN Jofek et al. 2007 : Smith et al. 2007 : Smith &: McCray 
2OO7I : ISmith et al.lboOSbl . boioi iMiller et al.lboiObh . SN 2006gy peaked at My ~ -22 mag, 
corresponding to a luminosity of L ~ 3 x 10^'^ ergs/ s, and it radiated of order i^rad > 2 x 10^^ 



ergs in total ( lOfek et al.l 120071 : iMiller et al.ll2010bh . In a normal Type II SN, a very small 
fraction (< 1%) of the available energy is radiated. The initial energy from the shock is 
converted into kinetic energy by adiabatic expansion, leaving little to power the luminosity. 
Thus, a supernova ejecting Mej ~ IOMq at velocities of fej ~ 4000 km/s contains enough 
kinetic energy, i^kin — 2 x 10^^ ergs to power the transient, but it is only radiated on time 
sc ales of order 10^ years as the SN remnant develops. The scenario we present was proposed 
by ISmith fc McCravl J2007h for SN 2006gy with further elaborations in ISmith et all J2008bh 



and Miller et al 



J2010bh ~ 



Very massive stars (> SOM©) are known to undergo impulsive events ejecting significant 
fractions of the stellar envelope. The rnost f amous example is the eruption of rj Carina in 



1837-1857 (see iDavidson fc Humphrevslll997f). but also includes some o f the so-called "su 



pernova impostors" (e.g. ISmith fc Owockil l2006l : iPastorello et al.l 120071 ). The number and 
temporal distributions of these mass ejections relative to the final supernova is unknown, but 
there are several cases where one seems to have occurred shortly before the final supernova 
(e.g., ISmith et al.ll2007l ). Pair instabiht y SNe can also produ ce a sequence of mass ejections 
just prior to the death of the star (see IWoosley et al.l 120071 ). Let us suppose that the pro- 
genitor of SDWFS-MT-l was such a massive star and had two such eruptions in its final 
phases, roughly 300 and 5 years before the final explosion, each o f which ejected of orde r 
^sheii — IOMq. For a typical ejection velocity of 200 km/s (e.g. ISmith &: Owockil |2006| ). 



these shells would lie roughly 13000 and 200 AU from the star when it exploded. 

The role of the first shell is t o conv ert the kinetic energy of the explosion into ra- 
diation following ISmith fc McCrayl (|2007| ). Roughly speaking, the collision between the 
ejected material and the first shell can extract fraction Msheii/ (Mej + Mgheii) of the ki- 
netic energy and convert it to radiation. The distance to the shell is tuned so that the 
shock crossing time is comparable to the Thomson photon diffusion time, thereby allow- 
ing the radiation to escape without significant adiabatic losses. This requires a total mass 
of approximately M = 47r/2^mpc/(jTfej = 2. 4000^0 where R = 100i?ioo AU and 
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fej = 4000fej,4ooo km/s. We can get up to 20Mq either by making R ~ 300 AU or by allowing 
for adiabatic losses. This then implies characteristic time, energy, luminosity and temper- 
ature scales of t = R/vej = 50i?ioo/^^cj,4ooo days, E = (l/4)Mfgj = 2 x 10^°-R^QQi;ej^4ooo ergs, 
E/t = 5x 10^3^2.^^^^^^^^ g^gg/g^ y ^ (L/47ri?V)i/^ = MOOOt^J^o^ror K, respec- 
tively, where we have used the estimate of the total mass and equally divided it between the 
shell and the SN ejecta. The b asic energetics, but no t the radiation transport, are largely 



confirmed by the simulations of Ivan Marie et al.l (120091 ) . This set our choice of incident tem- 
perature in the DUSTY models. Thus, a shell of ejected material at i? ~ SOOfgj 4QQg AU 
can produce the necessary energetics and luminosities but cannot match the observed time 
scales (too short) or temperatures (too high). 



The rol e of the second shell is to match the time and temperature scales. ISmith et al. 



fl2008tf l and Miller et all fl2010tf l invoke a second M^heii ^ IOMq shell located at ~ 1 pc 
in order to explain the relatively bright, long-lived IR emission observed for SN 2006gy. 
Their shell is, however, optically thin to optical light (r 0.02) because the total infrared 
emission is only a small fraction of the total. The optical depth of a shell is of order 
T ~ 0.1fi:(Msheii/10MQ)(10^AU/i?)^ for k in cm^/g, which leads to an optically thin shell with 
Topt — 0.08 for an R 1 pc shell radius a nd ^opt — 500 cm^/g for a radiation temperature 



of order 10^ K (e.g., ISemenov et al.l 120031 ). This opacity, n ~ 500 cm^/g, corresponds to 
a dust-to-gas ratio of approximately 1%, which arguably may be high for a low metallicity 
galaxy. Note, however, that SN 2006gy is probably also associated with a low metallicity 
galaxy (see §3.4) and the shell producing its dust echo contains enough dust for our scenario. 

To explain SDWFS-MT-1 we simply move the shell inward. At this point, however, we 
should improve on Eqn. ([1]) with a simple model for an optically thick dust echo. Suppose 
that the dust lies in an optically thick spherical shell of radius R and is exposed to luminosity 
Lp for a period of time tp. Because all the dust is exposed to the same luminosity, the dust 
temperature is independent of position or time while irradiated. The observer, however, sees 
the luminosity pulse modulated by the light travel time 2R/c across the shell. In this simple 
model, the light curves are trapezoids with linear rises and falls connected by a plateau 
created by the changing fraction of the shell illuminated by the outburst from the observer's 
perspective, but the dust temperature is constant. Models that can reproduce the drop 
in luminosity between epochs 3 and 4 have the pulse time tp < 2R/c, with the simplest 
model having the luminosity plateau run from epoch 2 to epoch 3. The plateau duration 
is 2R/c — tp = 180 days while the rise and falls last tp. The luminosity of the plateau 
is 27rRctpCrT^, which is smaller than Lp = AirR'^crT^ by the ratio 2tp/Rc ~ 0.2 between 
the light crossing time and the pulse length (if tp > 2R/c the plateau luminosity is Lp). 
The factor of 3 drop in luminosity between epochs 3 and 4 is L4/L3 = 1 — t34/tp ^ 1/3 
where ^34 ~ 33 days is the time between these epochs. Thus, we have tp ~ 50 days. 
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and R ~ 19000 AU or 115 light days. With these time and distance scales we come close to 
matching the epoch 2 and 3 luminosities given the observed temperatures, with an estimated 
luminosity of L2 = L3 = 5 x 10^^ ergs/s. The true luminosity is Lp ~ 2 x 10^^ ergs/s, which 
is why the radius estimate is somewhat larger than the estimates of Eqn. ([T]). There is some 
tension between the plateau luminosity and the drop in luminosity between the last two 
epochs, but the overall agreement is good given the crudeness of the model both for the dust 
emission and the luminosity profiles. 

The mass required to produce an a shell with optical depth r at this distance is M ~ 
r(K/500cm^/g)(i?/19000 AU^Mq, so there is little difficulty having an optically thick shell 
for masses of order I-IOM© even if the dust-to-gas ratio needs to be significantly reduced from 
~ 1% because of the low nietallicity of the galaxy . While the dust temperature is close to that 



for destroying dust (e.g., iDwek fc Wernerlll98ll ). the situation is somewhat different from a 



normal supernova because the (Thomson) optically thick inner shell protects the outer from 
any initial high luminosity bursts of hard radiation. For the IR-radiation, kjr ~ 10 cm^/g, 
and the shell is optically thin to the mid-IR emission for reasonable shell masses, so the dust 
emission streams freely to the observer. Given the estimated pulse duration, tp ~ 50 days, 
and luminosity, Lp = 2 x 10^^ ergs/s, we can estimate that the inner shell is at i? ~ 160 AU, 
the velocity must be Wej — 6300 km/s and the mass involved must be M ~ 4:Mq, where the 
rapid decline pushes these values to small radii/masses and high velocities in order to keep 
tp short. Given the crudeness of the model, it is remarkable how well all the observational 
features can be matched using relatively sensible physical parameters. 



3.3. The Transient's Host Galaxy 



The nature of the host galaxy indirectly supports this hypothesis because both GRBs 
and most hyper-luminous Type Iln SNe are also assoc iated with low lumi nosity, low metal- 
licity dwarf galaxies. Figure [5l similar to Figure 1 in IStanek et al.l (120061 ) . shows the posi- 
tion of the host galaxy for SDWFS-MT-1 i n the metallicity-luminosity plane in relation to 



126 , 000 star-forming S PSS DR4 galaxies (ITremonti et al.ll2004j : lAdelman-McCarthy et al. 



20061 : iPrieto et al.l 120081 ). We also include the hosts of six local {z < 0.25) GRBs using 



updated estimates from iLevesaue et al.l (120101) and including the host of GRB 100316D 
Jchornock et al.lboiol : Istarling et allboioh . While SDWFS-MT-1 was not a GRB, its oc- 
currence in a very low metallicity galaxy is striking, since such galax ies produce only a s mall 
fraction of massive stars in the local universe (e.g., see Figure 2 in IStanek et al.ll2006l ). In 
Figure O we ignore small (< 20% ) effects on metallicity measurements arising from the 
method being used to derive it (see lKewley &: Ellisonll2008l ). 
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In fact, the hosts of most other hyper-luminous SNe also tend to be very low luminosity 
(~ 0.01 L^), suggestive of their low metallicity (unfortur iately, oxygen abund ances have not 
been measured for most of these hosts). SN 2005ap (iQuimby et al.l 120071 ) a.t z = 0.283 
peaked at unfiltered magnitude of M = —22.7 mag and is associated with a Mr = —1 6.8 
{Mb ^ -16.1) mag dwarf galaxy. SN 2006tf JOuimbv et al.ll2007l : ISmith et al.[l2008ah at 
z = 0.074 exploded in a host of low luminosity Mr = —16.9 {Mb ~ —16.2) mag. The 
over- luminous M/j = —21.3 mag Type Ic SN 200 7bi exploded in a fain t , Mr = —16 . 3 ma g 
dwarf galaxy, with 12 + l og(0/H) = 8.18 ± 0.17 JCal-Yam et!J l2009|; lYoung et al l boioh . 
SN 2008es at z = 0.213 f lMiller et allbood : ICezari et aLlboOoh reached My = -22.3 mag 
and exploded in a host which was not detected down to a faint luminosity of My > —17.4 
{Mb > -17.1) mag. SN 2008fz at z = 0.133 peaked at My = -22.3 m ag, but no host ha s 
been detect ed down to an upper limit of Mr > —17 (Mn > — 16.6) mag (jPrake et al.ll2009l ). 
SN 2006gy jQuimbvl[2006l : ISmith et al. l l2007l : lofek et ahlEoOTl ) is one of two exceptions here, 
as it exploded in a fairl y luminous SO g alaxy, NGC 1260 with Mb ~ —20.3 mag. Another 
example is SN 2003ma (IRest et al.l 120091 ) that released 4 x 10^^ ergs in a Mb ~ —19.9 mag, 
12 + log(0/H) ^ 8.7 galaxy. 

The general message is clear: GRBs and luminous SNe prefer low metallicity, low lumi- 
nosity galaxies, so does SDWFS-MT-1. 



3.4. Rates 

We systematically searched for transients similar to SDWFS-MT-1 and SNe in general 
in the SDWFS data. We considered highly correlated [3.6] and [4.5] light curves (r > 0.8) 
with flag = 1 (photometry unaffected by the wings of bright stars), variabi lity significances 



crp.fi] > 1 and cr[4.5] > 1 and colors [3.6] — [4.5] < 0.4 to eliminate AGN (see iKozlowski et al. 



2010al ). We selected as candidates objects that at their peak were brighter than [3.6] < 
18.67 mag, corresponding to 1 mag brighter than the 3a detection limit, and which changed 
in brightness by at least one magnitude between either epochs 1 and 2 or epochs 2 and 3/4. 
The latter two epochs are close enough together that an SN might not show a significant 
brightness change between them. We found six candidate objects, one of which is SDWFS- 
MT-1. The next four of them are artifacts due to PSF structures (long spikes) from bright 
stars. 

The remaining object, SDWFS J142557.64+330732.6 (hereafter SDWFS-MT-2), peaks 
in epoch 2 close to our selection limit at [3.6] = 18.52±0.14 mag and [4.5] = 18.71±0.33 mag. 
It lies 0'.'7 South- West of a potential host galaxy with GALEX NUV = 24.02 ± 0.19 mag, 
NDWFS 5w = 23.77 ± 0.03 mag, R = 23.34 ± 0.04 mag, and / = 22.68 ± 0.06 mag that is 
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not detected in SDWFS. In the deep images from NDWFS, the host galaxy appears to be an 
extended, possibly edge-on (length 3.0 arcsec in 1.2 arcsec seeing) disk, with a redder core and 
bluer outer regions. The transient lies in the bluer outer regions, although its position is only 
known to ~ 0'.'2 due to its faintne ss. We estimate a photometric redshift of 2; 1.0 for the 
host using the template models of lAssef et al.l ( l2010l ). corresponding to a comoving distance 
of 6.6 Gpc and a linear scale of 8.0 kpc/arcsec. This translates into the linear diameter of 
the host of ~ 24 kpc. The best fit model for the SED of the host is a mixture of an AGN 
and a starburst template. The mid-IR color of the transient, [3.6] — [4.5] = 0.2 ± 0.5 mag, 
is ambiguous on the question of an SN or AGN origin for the transient, although the offset 
position would favor an SN. Given the redshift estimate, however, the impli ed luminosity is 
too h igh even for the brightest normal SN - for M[3 g] ~ Mk ~ —20 mag (IMannucci et al. 
2OO3I ) the peak magnitude of [3.6] ~ 24 mag would be far below our detection limits. A 



second event like SWFS-MT-1 would have to be even more extreme. Given the paucity of 
the data, interpreting this source as a false positive seems the most probable explanation. 

For events with six month durations like SDWFS-MT-1 our survey time coverage is 
approximately AT = 2 years. SDWFS-MT-1 peaked in epoch two at fluxes of 0.12 mJy in 
[3.6] and 0.10 mJy in [4.5] at a luminosity distance of 920 Mpc. The peak flux limit for 
our systematic search is 9.6/iJy (5(T in [4.5]) which means we could have detected SDWFS- 
MT-1 at a luminosity distance of 3000 Mpc (920 x (100/9.6)^/^) about z = 0.52 ignoring 
/^-corrections, which will initially aid detection because the SED peaks bluewards of the 
IRAC bands. Thus, we could detect the source in a comoving volume of approximately 
V = 6 X 10^ Mpc^. Combining these factors, we estimate a rate for SDWFS-MT-1- like 
transients of r = {VAT)~^ ~ 8 x 10~^ Mpc^^ yr~^ where we have not corrected the time 
scales for redshifting effects given the otherwise large uncertain ties from the statistic s of one 
event. This is approximately 10~^ of the Type II SN e rate (IHoriuchi et al.l 120091 ), ~ 0.1 
of th e hyper-luminous Type II SNe (IMiller et al.ll2009l ). and ~ 10 of the GRB rate (INakar 
2007h . 



4. Summary 

In this paper, we presented the discovery of a luminous mid-IR transient, SDWFS- 
MT-1, that is unusual in several ways. Its most striking features are the high luminosity, 
L ~ 6 X 10'^^ ergs/s, the low apparent temperature, T ~ 1350 K, the long, 6 month duration 
of the event, and the resulting high amount of radiated energy, Ej-^d ^ 10^^ ergs. A model 
that explains the transient is a self-obscured supernova, that has two concentric (patchy) 
shells, each with the mass of ~ 10 Mq, ejected prior to the SN explosion. The inner shell 
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at ~ 150 AU protects the dust in the outer shell from initial luminosity spikes and converts 
the high kinetic energy of the SN shock wave into thermalized radiation with T 10^ K. 
The dust in the more distant (0.05 pc) second shell absorbs the radiation of the first shell, 
and reradiates it at T f» 10^ K, while simultaneously smoothing out and stretching the light 
curve. This scenario matches the observed energy, luminosity, time and temperature scales 
of the tra nsient by slightly rn o difying the scenarios used for other hyper-luminous Type Iln 
SNe (e.g.. lSmith et al.ll2008al Jbl: iMiller et al.ll2010bl ). The host galaxy of the transient is also 
typical of these SNe, and our estimate of the rate is consistent with such events being ~ 10% 
of hyper-luminous Type II SNe. 



The discovery of the progenitors of SN 2008S and NGC300-OT-1 flPrieto et al.l 12008 
Thompson et al.ll2009l ) as self-obscured (T ~ 400 K, L ~ IO^Lq) stars already indicated that 
mid-IR observations are necessary to understanding the fates of massive stars. If our inter- 
pretation of SDWFS-MT-1 is correct, then this may extend to the actual transients as well as 
to their progenitors, although it is quite likely that SDWFS-MT-1 w as also a significant op 



tical transient. IR dust e chos from SNe are relatively common, (e.g.. lDweklll983l : iRest et al. 



20091 : iMiller et al.ll2010bl ). but in these cases the optical depth of the dust is low and the IR 
emission is only a fraction of what is o bserved in the optica l. A few recen t SNe with signifi- 



cant echo luminosit ies are SN 20 07od ([Andrews et al. 2010 ). SN 2006gy ( Smith et a. 



Miller et al.l l2010b[ ). SN 2004et (IKotak et al.ll2009h . and SN 2002hh (iMeikle et al 



2007 



20061 ). 



Obviously there is a strong selection effect against optically selecting SNe where the opti- 
cal depth of the surviving dust is near unity, so the true abundance of such sources can 
only be deter mined in near/raid-IR searches. There h ave been some searches made in the 
near-IR (e.g., iMannucci et al.ll2003l : ICresci et al.l 120071 ). principally focused on the problem 
of unassociated foreground absorption rather than dust associated with the progenitor star. 

Two important steps are to try and identify massive sta rs with such deri s e she lls of 
ejected mat e rial, s imilar to the search for self-obscured stars in iThompson et al.l ( l2009l ) and 
Khan et al.l (120101 ). and to better characterize the mid-IR emission of other hyper-luminous 
Type Iln SNe. If the inner shell invoked to convert the kinetic energy into radiation is dusty, 
the progenitor will appear as a rapidly cooling near/mid-IR source as the shell expands. 
Such shells will have temperatures of order the 400 K observed for SN 2008S or NGC300- 
OT- 1 but should show the time variability known to be absent for the SN 2008S progenitor 
(see iPrieto et al.ll2008l ). For the dust to survive the SN explosion and produce a powerful 
mid-IR transient, it must be considerably more distant and colder than the T ~ 400 K dust 
photosphere of the progenitors of SN 2008S or NGC300-OT-1, with temperatures of order 
100 K and SEDs peaking at 20-30/im that will make them invisible to warm Spitzer. In fact, 
the progenitors may well resemble a moderately colder version of rj Carina, which radiates 
most of its energy in the mid-IR and has an SED peaking at ~ 10/im (jRobinson et al. 
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19871 ). The temperature of such a distant shell will evolve slowly because the time scale for 
the expansion of the shell is now quite long. Warm Spitzer will remain a valuable tool for 
identifying warm, inner shells, but a careful search for cold outer shells will need to wait for 
the James Webb Space Telescope (JWST). 



The Wide-Field Infrared Survey Explorer ( WISE, iMainzer et al.l l2005l : I Wright et al. 



20101), launched in December 2009, is currently carrying out an all-sky survey at similar 
wavelengths to Spitzer, as a natural by-product of its design, find many transients sim- 
ilar to SDWFS-MT-1. With a 5a se nsitivity of 0.12 mJ y at 3.4 fim and 0.16 mJy at 
4.7 fim for a region with 8 exposures ( IMainzer et al.ll2005l ). WISE has far less sensitivity 
than the SDWFS survey. Crudely speaking the survey volume compared to SDWFS scales 
as the ratio of the {solid angle) {sensitivity)~^^'^ of the two surveys, so Vwise/Vsdwfs — 
(42000/8)(0.12mJy/6.3;uJy)~3/2 _ ^q2 greater sohd angle wins over the reduced sen- 

sitivity. If WISE completes two full years of observations it will also have 4 epochs, 1 epoch 
every 6 months, and so should detect a significant number of these mid-IR transients (~ 10^) 
independent of their origin. For sources closer to the poles of the WISE scan pattern, WISE 
will be able to produce detailed hght curves. On longer time scales, JWST will be able to 
identify and monitor such luminous transients at almost any redshift. 
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Epoch 1 Epoch 2 Epoch 3 Epoch 4 

• 4* ca i| 




Fig. 1. — SDWFS images, approximately 1.5 x 1 arcmin^, centered on the transient. The rows 
correspond to the [3.6], [4.5], [5.8], and [8.0] bands (downward from top), and the columns 
show Epochs 1 through 4 (left to right), taken in 2004, 2007 and twice in 2008, respectively. 
The object is marked with the cross hairs. The source is clearly absent in Epoch 1 and then 
present in all subsequent epochs except for the Epoch 4 [8.0] data. 
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Fig. 2. — Light curves of the transient. Top paneh observed optical and mid-IR magnitudes. 
For display purposes the IRAC points are offset by —3 days [3.6], —1 day [4.5], +1 day [5.8] 
and, +3 days [8.0]. The right y-axis shows the absolute magnitudes without fT-corrections. 
A normal Type II SN would peak at M ~ —20 mag at 3.6;um. The gray area marks the 
region when the mid-IR transient was not visible to airmass < 1.5. The large dark green 
symbols mark the CRTS upper limits found by stacking the three nearest epochs (small 
symbols) into one image. The RAPTOR upper limits (purple symbols) are stacks of 21-62 
images. The QUEST upper limits (orange) are derived based on 1-4 stacked images per 
epoch. The large symbol for QUEST combines the four nearest epochs. Bottom panel: 
points from the upper panel converted to luminosity as \Fx. 
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Fig. 3. — The Kcck/LRIS rest-frame spectrum (selected parts) of the transient's host galaxy 
with the major emission lines marked. From the line positions we find a redshift of z = 
0.1907. To derive the metallicity, we measured the fluxes of the [Oil] 3727A, [OIII] 4363A, 
[OIII] 4959A, [OIII] 5007A and H/3 486lA lines. 
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Fig. 4. — SED of the transient during the second epoch of SDWFS. We fit a black body SED 
that includes the dust emissivity {Q oc A~^) to the four IRAC bands (filled circles) and find a 
best fitting temperature of T = 1357 ±57 K (dotted l ine). We also sho w a template SED for 
the host galaxy modeled as an irregular (dashed line; Assef et "aPboiol ) and the combination 
of the modified black body and irregular galaxy SEDs (thick solid line). Overplotted are 
the V, R and IRAC filter transmission curves, "normalized" to 5 x 10^^ ergs/s. The MDM 
(QUEST) limit is derived for an observation taken 90 (~ 40) days prior to the second epoch 
of SDWFS. The CRTS limit corresponds to the stack of 3 epochs (12 x 30 sec images) taken 
between SDWFS epochs 3 and 4. The open squares mark the NDWFS Ri and /-band 
measurements for the host galaxy (from left to right). 
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Fig. 5. — Metallicity-Iuminosity relatio n for 125,958 SDSS DR4 star forming galaxies (e.g, 



Tremonti et al.ll2004l : iPrieto et al.l 120081 ) . They are binned into 0.1 mag and 0.025 dex bins 
and the smoothed contours (3x3 bins) are drawn for 1, 10, 50, 100 and 20 objects per 



bin, co unting from the outer contour. The figure is analogous to Figure 1 from lStanek et al. 



(]2006f ). SDWFS-MT-1 is shown as the filled pentagon (lower left). A handful of local GRBs 
with associated SNe (and two SNe) with known hosts metallicities and luminosities are 
shown with open circles (small pentagons). The arrows show Mb for five host galaxies of 
high luminosity SNe with unknown oxygen abundances. They are, from left to right, SN 
2005ap, SN 2006tf, SN 2008fz, SN 2008es, and SN 2006gy (see Section B- The LMC/SMC 
oxyg en abundances are ado pted following iPeeples et al.l (|2008[ ) and Milky Way's (s o lar) a re 
from Delahaye et al. ( 2010 ). while the absolute magnitudes are from iKarachentsevI (2005). 
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Table 1. SDWFS Mid-IR Data In The Observed Frame. 



Band (Units) 


Epoch 1 


Epoch 2 


Epoch 3 


Epoch 4 


\3.6] fmag) 


> 19.67 


15.93 ±0.01 


16.09 ±0.01 


17.99 ± 0.08 


[4.5] (mag) 


> 18.73 


15.61 ±0.02 


15.76 ±0.03 


17.72 ±0.15 


[5.8] (mag) 


> 16.33 


15.44 ±0.15 


15.74 ±0.21 


16.93 ±0.55 


[8.0] (mag) 


> 15.67 


15.43 ±0.30 


15.09 ±0.26 


> 15.67 


AFa([3.6]) (10^0 ergs/s) 


< 32 


1017 ±18 


878 ± 15 


153 ± 11 


AFa([4.5]) (10^° ergs/s) 


< 39 


690 ± 16 


601 ± 19 


99± 14 


AFa([5.8]) (10^° ergs/s) 


< 179 


406 ± 42 


308 ± 61 


103 ± 55 


AFx([8.0]) (10^° ergs/s) 


< 133 


166 ± 47 


227 ± 55 


< 133 




Black body fit including 


dust emissiv 


ty {Q oc A"^) 






^rad > 0.9 X 10^1 ergs 




Black Body Temp. (K) 




1357 ± 57 


1338 ± 71 


1400 ± 353 


"Total" L {10^^ ergs/s) 




~ 6.8 


~ 5.7 


~ 1.8 


Radius R (AU)« 




~ 11000 


~ 10500 


~ 5500 



Note. — The SDWFS magnitudes were measured in 4 arcsec diameter aperture 
with aperture corrections to 24 arcsec. 

"Radius R corresponds to the radius where the optical depth Topticai ~ 1- 



